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Abstract 

Consider an infinite collection of qubits arranged in a line, such that every 
pair of nearest neighbors is entangled: an "entangled chain." In this paper 
we consider entangled chains with translational invariance and ask how large 
one can make the nearest neighbor entanglement. We find that it is possible 
to achieve an entanglement of formation equal to 0.285 ebits between each 
pair of nearest neighbors, and that this is the best one can do under certain 
assumptions about the state of the chain. 



PACS numbers: 03.67.-a, 03.65.Bz, 89.70. +c 



1 Introduction: Example of an entangled chain 



Quantum entanglement has been studied for decades, first because of its im- 
portance in the foundations of quantum mechanics, and more recently for its 
potential technological applications as exemplified by a quantum computer 
[1]. The new focus has led to a quantitative theory of entanglement [2, 3, 4] 
that, among other things, allows us to express analytically the degree of en- 
tanglement between simple systems [5]. This development makes it possible 
to pose new quantitative questions about entanglement that could not have 
been raised before and that promise fresh perspectives on this remarkable 
phenomenon. In this paper 1 would like to raise and partially answer such 
a question, concerning the extent to which a collection of binary quantum 
objects (qubits) can be linked to each other by entanglement. 

Imagine an infinite string of qubits, such as two-level atoms or the spins 
of spin- 1/2 particles. Let us label the locations of the qubits with an inte- 
ger j that runs from negative infinity to positive infinity. I wish to consider 
special states of the string, satisfying the following two conditions: (i) each 
qubit is entangled with its nearest neighbors; (ii) the state is invariant under 
all translations, that is, under transformations that shift each qubit from its 
original position j to position j + n for some integer n. Let us call a string of 
qubits satisfying the first condition an entangled chain, and if it also satisfies 
the second condition, a uniform entangled chain. Note that each qubit need 
not be entangled with any qubits other than its two nearest neighbors. In this 
respect an entangled chain is like an ordinary chain, whose links are directly 
connected only to two neighboring links. By virtue of the translational in- 
variance, the degree of entanglement between nearest neighbors in a uniform 
entangled chain must be constant throughout the chain. The main question 
I wish to pose is this: How large can the nearest-neighbor entanglement be 
in a uniform entangled chain? 

This problem belongs to a more general line of inquiry about how entan- 
glement can be shared among more than two objects. Some work on this 
subject has been done in the context of the cloning of entanglement [6—11], 
where one finds limits on the extent to which entanglement can be copied. In 
a different setting not particularly involving cloning, one finds an inequality 
bounding the amount of entanglement that a single qubit can have with each 
of two other qubits [12] . One can imagine more general "laws of entanglement 
sharing" that apply to a broad range of configurations of quantum objects. 
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The present work provides further data that might be used to discover and 
formulate such laws. The specific problem addressed in this paper could also 
prove relevant for analyzing models of quantum computers in which qubits 
are arranged along a line, as in an ion trap [13]. The infinite chain can be 
thought of as an idealization of such a computer. Moreover, the analysis of 
our question turns out to be interesting in its own right, being related, as we 
will see, to a familiar problem in many-body physics. 

To make the question precise we need a measure of entanglement between 
two qubits. We will use a reasonably simple and well-justified measure called 
the "concurrence," which is defined as follows [5]. 

Consider first the case of pure states. A general pure state of two qubits 
can be written as 

1^) =«|00)+/?|01)+7|10) + 5|11). (1) 

One can verify that such a state is factorizable into single-qubit states — that 
is, it is unentanglcd — if and only if q5 = f3j. The quantity C = 2\a6 — f3'~f\, 
which ranges from to 1, is thus a plausible measure of the degree of entan- 
glement. We take this expression as the definition of concurrence for a pure 
state of two qubits. For mixed states, we define the concurrence to be the 
greatest convex function on the set of density matrices that gives the correct 
values for pure states [14]. 

Though this statement defines concurrence, it does not tell us how to 
compute it for mixed states. Remarkably, there exists an explicit formula for 
the concurrence of an arbitrary mixed state of two qubits [5]: Let p be the 
density matrix of the mixed state, which we imagine expressed in the stan- 
dard basis {|00), |01), |10), |11)}. Let p, the "spin-flipped" density matrix, be 
{ay ® ay)p*{ay (g) ay), where the asterisk denotes complex conjugation in the 

standard basis and ay is the matrix ^ ^ 0^ ) Finally, let Ai, A2, A3, A4 be 

the square roots of the eigenvalues of pp in descending order — one can show 
that these eigenvalues are all real and non-negative. Then the concurrence 
of p is given by the formula 

C{p) = max{Ai - A2 - A3 - A4, 0}. (2) 

The best justification for using concurrence as a measure of entanglement 
comes from a theorem [5] showing that concurrence is a monotonically in- 
creasing function of the "entanglement of formation," which quantifies the 
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non-local resources needed to create the given state [4].^ As mentioned above, 
the values of C range from zero to one: an unentangled state has C = 0, 
and a completely entangled state such as the singlet state ;^(|01) — |10)) 
has C — 1. Our problem is to find the greatest possible nearest-neighbor 
concurrence of a uniform entangled chain. At the end of the calculation we 
can easily re-express our results in terms of entanglement of formation. 

Another issue that needs to be addressed in formulating our question is 
the meaning of the word "state" as applied to an infinite string of qubits; in 
particular we need to discuss how such a state is to be normalized. Formally, 
we can define a state of our system as follows. A state w of the infinite string 
is a function that assigns to every finite set S of integers a normalized {i.e., 
trace one) density matrix w{S), which we interpret to be the density matrix 
of the qubits specified by the set S; moreover the function w must be such 
that if 5*2 is a subset of 5*1, then 10(82) is obtained from w{Si) by tracing 
over the qubits whose labels are not in S2. This formal definition is perfectly 
sensible but somewhat bulky in practice. In what follows we will usually 
specify states of the string more informally when it is clear from the informal 
specification how to generate the density matrix of any finite subset of the 
string. We will also usually use the symbol p instead of w{S) to denote the 
density matrix of a pair of nearest neighbors. 

It is not immediately obvious that there exists even a single example of 
an entangled chain. Note, for example, that the limit of a Schrodinger cat 
state — an equal superposition of an infinite string of zeros with an infinite 
string of ones — is not an entangled chain. In the cat state, the reduced 
density matrix of a pair of neighboring qubits is an incoherent mixture of 
|00) and |11), which exhibits a classical correlation but no entanglement. 
(Note, by the way, that our informal statement "an equal superposition of 
an infinite string of zeros with an infinite string of ones," specifies exactly 

^One can define the entanglement of formation as follows. Let p be a mixed state 
of a pair of quantum objects, to be shared between two separated observers who can 
communicate with each other only via classical signals. The entanglement of formation 
of p is the asymptotic number of singlet states the observers need, per pair, in order to 
create a large number of pairs in pure states whose average density matrix is p. (This 
is conceptually different from the regularized entanglement of formation, which measures 
the cost of creating many copies of the mixed state p [15]. However, it is conceivable that 
the two quantities are identical.) Entanglement of formation is conventionally measured 
in "ebits," and for a pair of binary quantum objects it takes values ranging from to 1 
ebit. 
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the same state as if we had taken an incoherent mixture of these two infinite 
strings: no finite set of qubits contains information about the phase of the 
superposition.) 

We can, however, construct a simple example of an entangled chain in 
the following way. Let Wq be the state such that for each even integer j, the 
qubits at sites j and j + 1 are entangled with each other in a singlet state. 
We can write this state informally as^ 

^/|0)_2|l)-l - |l)-2|0)_l\ /|0)o|l)l - |l)o|0)i\ /|0)2|1)3 - |1)2|0)3\ 



V V2 J \ J \ V2 J 

(3) 

The state Wq is not an entangled chain because the qubits are not entangled 
with both of their nearest neighbors: qubits at even-numbered locations are 
not entangled with their neighbors on the left. However, if we let Wi be the 
state obtained by translating wq one unit to the left (or to the right — the 
result is the same), and let w be an equal mixture of Wq and Wi — that is, 
w — {wq + Wi)/2 — then w is a uniform entangled chain, as we now show. 

That w is translationally invariant follows from the fact that both wq and 
Wi are invariant under even displacements and that they transform into each 
other under odd displacements. Thus we need only show that neighboring 
states are entangled. For definiteness let us consider the qubits in locations 
j — 1 and j — 2. In the state wq, the density matrix for these two qubits is 



that is, the completely mixed state. (The two qubits are from distinct singlet 

^Alternatively, we can characterize the state wq according to our formal definition 
by specifying the density matrix of each finite collection of qubits: Let S define such 
a collection. Then for each even integer j such that both j and j + 1 are in S, the 
corresponding pair of qubits is in the singlet state; all other qubits {i.e., the unpaired 

f - \ 

ones) are in the completely mixed state |^ ^ i J ' ^^'^ ^^^^ density matrix w(<S') is 
obtained by taking the tensor product of the pair states and single-qubit states. 
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pairs.) The density matrix of the same two qubits in the state wi is 



/ 

° K 

V 



0\ 



0/ 



1 

2 





(5) 



that is, the singlet state. In the state w, the qubits are in an equal mixture 
of these two density matrices, which is 



I 



\ 




V U 



-i f 



(6) 



It is easy to compute the concurrence of this density matrix, because p is the 
same as p itself. The values Aj in this case are the eigenvalues of p, which are 
|, |, |, |. The concurrence is therefore C— ^ — \ — \ — \ — \. This same 
value of the concurrence applies to any other pair of neighboring qubits in the 
string because of the translational invariance. The fact that the concurrence 
is non-zero implies that neighboring qubits are entangled, so that the state 
w is indeed an entangled chain. For uniform entangled chains, we will call 
the common value of C for neighboring qubits the concurrence of the chain. 
Thus in the above example the concurrence of the chain is |. 

As we will see, it is possible to find uniform entangled chains with greater 
concurrence. Let Cmax be the least upper bound on the concurrences of all 
uniform entangled chains. We would hke to find this number. We know that 
Cmax is no larger than 1, since concurrence never exceeds 1. In fact we can 
quickly get a somewhat better upper bound, using the following fact: when 
a qubit is entangled with each of two other qubits, the sum of the squares of 
the two concurrences is less than or equal to one [12]. In a uniform entangled 
chain, each qubit must be equally entangled with its two nearest neighbors; 
so the concurrence with each of them cannot exceed l/v^. Thus, so far what 



we know about Cmax is this: 



1/4 < Cmax < 1/V2. 



(7) 



This is still a wide range. Most of the rest of this paper is devoted to getting 
a better fix on Cmax by explicitly constructing entangled chains. 
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2 Building chains out of blocks 



Using the above example as a model, we will use the following construction 
to generate other uniform entangled chains. (1) Break the string into blocks 
of n qubits, and define a state wq in which each block is in the same n-qubit 
state |(^); that is, wq is a tensor product of an infinite number of copies of |^). 
(In the above example n had the value 2 and |^) was the singlet state.) (2) 
Define w^, k — 1, . . . ,n — 1, to be the state obtained by shifting wq to the 
left by k units. (3) Let the final state w be the average {wq + ■ — h Wn-i)/n. 
A state generated in this way will automatically be translationally invariant. 
In order that the chain have a large concurrence, we will need to choose the 
state 1^) carefully. Finding an optimal |^) and proving that it is optimal may 
turn out to be a difficult problem. In this paper I will choose |^) according 
to a strategy that makes sense and may well be optimal but is not proven to 
be so. 

In the final state w, each pair of neighboring qubits has the same density 
matrix because of the translational invariance. Our basic strategy for choos- 
ing 1^), described below, is designed to give this neighboring-pair density 
matrix the following form: 



P 



( pll 













P22 


P23 








P23 


P33 





I 








/ 



(8) 



(The ordering of the four basis states is the one given above: 1 00) , 1 01) , 1 10) , 1 11) .) 
One can show that the concurrence of such a density matrix is simply 

C = 2|p23|. (9) 

Besides making the concurrence easy to compute, the form (8) seems a rea- 
sonable goal because it picks out a specific kind of entanglement, namely, a 
coherent superposition of |01) and |10), and limits the ways in which this en- 
tanglement can be contaminated or diluted by being mixed with other states. 
In particular, the form (8) does not allow contamination by an orthogonal 
entangled state of the form q;|00) + — orthogonal entangled states when 
mixed together tend to cancel each other's entanglement — or by the combi- 
nation of the two unentangled states |00) and |11). If the component P44 were 
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not equal to zero and the form were otherwise unchanged, the concurrence 
would be C = max{2(|p23| — y^pTTpil), 0}; so it is good to make either pu or 
P44 equal to zero if this can be done without significantly reducing p23- We 
have chosen to make P44 equal to zero. 

As it happens, one can guarantee the form (8) for the density matrix of 
neighboring qubits by imposing the following three conditions on the n-qubit 
state 1^): (i) |{) is an eigenstate of the operator that counts the number of 
qubits in the state |1). That is, each basis state represented in |^) must 
have the same number p of qubits in the state |1). (ii) \^) has no component 
in which two neighboring qubits are both in the state |1). (iii) The nth 
qubit is in the state |0). (This last condition effectively extends condition 
(ii) to the boundary between successive blocks.) Condition (i) guarantees 
that the density matrix p for a pair of nearest neighbors is block diagonal, 
each block corresponding to a fixed number of I's in the pair. That is, there 
are two single-element blocks corresponding to |00) and |11), and a 2x2 block 
corresponding to |01) and |10). Conditions (ii) and (iii) guarantee that P44 
is zero. The conditions thus give us the form (8). We impose these three 
conditions because they seem likely to give the best results; we do not prove 
that they are optimal. 

To illustrate the three conditions and how they can be used, let us consider 
in detail the case where the block size n is 5 and the number p of I's in each 
block is 2. (Our strategy does not specify the value of either n or p; these 
values will ultimately have to be determined by explicit maximization.) In 
this case, the only basis states our conditions allow in the construction of |^) 
are |10100), |10010), and |01010). Any other basis state either would have a 
different number of 1 's or would violate one of conditions (ii) and (iii) . Thus 
we write 

\C) -= aiallOlOO) + ai4|10010) + a24|01010). (10) 

The subscripts in aij indicate which qubits are in the state |1). The state w 
of the infinite string is derived from |0 as described above. We now want 
to use Eq. (10) to write the density matrix p of a pair of nearest neighbors 
when the infinite string is in the state w. For definiteness let us take the two 
qubits of interest to be in locations j = I and j = 2, and let us take the 
5-qubit blocks in the state Wq to be given by j = 1, . . . , 5, j = 6, . . . , 10, and 
so on. Our final density matrix p will be an equal mixture of five density 
matrices, corresponding to the five different displacements of Wq (including 
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the null displacement). 

For wq itself, the qubits at j = 1 and j = 2 are the first two qubits of |^). 
The density matrix for these two qubits, obtained by tracing out the other 
three qubits of the block, is 



P 



(0) 
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|ai3p + |ai4p 
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0/ 



(11) 



For wi, the qubits at j = 1 and j — 2 are now the second and third qubits 
of the block, since the block has been shifted to the left. Thus we trace over 
the first, fourth, and fifth qubits to obtain 
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(12) 



In a similar way one can find p^^^ and p^^^: 
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The density matrix corresponding to ^4 is different in that the two relevant 
qubits now come from different blocks: the qubit at j = 1 is the last qubit 
of one block and the qubit j — 2 is the first qubit of the next block. The 
corresponding density matrix is thus the tensor product of two single-qubit 
states: 



.(4) 
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To get the neighboring-pair density matrix corresponding to our final state 
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we average the above five density matrices, with the following simple result: 
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(13) 



where 

X = a\^au + 014024- (14) 
According to Eq. (9), the concurrence of the pair is 

2 

— I (3^3(314 + (3^4(3241 • (1^) 


Continuing with this example — n — 5 and p — 2 — let us find out what 
values we should choose for ai3, ai4, and 024 in order to maximize C . First, 
it is clear that we cannot go wrong by taking each aij to be real and non- 
negative — any complex phases could only reduce the absolute value in Eq. (15) — 
so let us restrict our attention to such values. To take into account the nor- 
malization condition, we use a Lagrange multiplier 7/2 and extremize the 
quantity 

013014 + O14O24 - (7/2) (oi3 + al^ + O24). (16) 

Differentiating, we arrive at three linear equations expressed by the matrix 
equation 

/ 1 \ / ai3 \ / ai3 \ 

1 1 oi4 = 7 «14 • (17) 

V 1 / V O24 / V «24 / 

Of the three eigenvalues, only one allows an eigenvector with non-negative 
components, namely, 7 = -\/2. The normalized eigenvector is 



/ ai3 \ 
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ai4 
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V2 




\ 024 / 
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/ 



which gives C — -\/2/5 = 0.283. This is greater than the value 0.25 that we 
obtained in our earlier example. 
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Before generalizing this calculation to arbitrary values of n and p, we 
adopt some terminology that will simplify the discussion. Let us think of 
the qubits as "sites," and let us call the two states of each qubit "occupied" 
(|1)) and "unoccupied" (|0)). The states \^) that we are considering have a 
fixed number p of occupied sites in a string of n sites; so we can regard the 
system as a collection of p "particles" in a one-dimensional lattice of length 
n. Condition (ii) requires that two particles never be in adjacent sites; it 
is as if each particle is an extended object, taking up two lattice sites, and 
two particles cannot overlap. Thus the number of particles is limited by the 
inequality 2p <n. 



3 Generalization to blocks of arbitrary size 

We now turn to the calculation of the optimal concurrence for general n and 
p assuming our conditions are satisfied. It will turn out that this calculation 
can be done exactly. 

For any values of n and p, the most general form of |^) consistent with 
condition (i) is 



10= E 



a 



(19) 



:;i<--<jp 



where . . . , jp) is the state of n sites j = 1, . . . , n in which sites ji, . . . , jp 
are occupied and the rest arc unoccupied. Because of conditions (ii) and (iii), 
aji,...jp must be zero if two of the indices differ by 1 or if jp has the value n. 
The coefficients in Eq. (19) satisfy the normalization condition 

E i«.w.r = i- (20) 
ji<---<jp 

Going through the same steps as in the above example, we find that in the 
state w the density matrix of any pair of neighboring sites is 



P 



1 



n 



( n-2p 






\ 

p y Q 

y* p 

y 



where 



E E E 

q=l ji<-<jp j[<-<j^ 



Jl, 



■,3'p^y^'U^^ n ^jUir 



(21) 



(22) 
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Here 6 is the Kronecker delta, and we define to be zero if any two 

of the indices are equal. In words, y is constructed as follows: Let two 
coefficients ciji,...,jp and cij[,...,j^ be called adjacent if they differ in only one 
index and if the difference in that index is exactly one; then y is the sum of 
all products of adjacent pairs of coefficients, the coefficient with the smaller 
value of the special index being complex conjugated in each case. In the 
above example there were only two such products, 0^3014 and 0^4024; hence 
the form of Eq. (14). 

As before, the concurrence of the chain is equal to 2|p23|; that is, C = 
{2/n)\y\. We want to maximize the concurrence over all possible values 
of the coefficients that are consistent with conditions (ii) and (iii). These 
conditions are somewhat awkward to enforce directly: one has to make sure 
that certain of the coefficients aji,...jp are zero. However, this problem is 
easily circumvented by defining a new set of indices. Let fci = ji, k2 = J2 — 1, 
^3 = Ja ~ 2, and so on up to kp = jp — (p — 1), and let = The 

constraints on the new indices kr are simply that < /ci < A:2 < ■ ■ ■ < A;^ < n', 
where n' = n — {p — 1). Finally, in place of |^), define a new vector \0: 

10= E h,,...,kp\h,...,K), (23) 

ki<---<kp 

where |A;i, . . . , kp) is the state of a lattice of length n' — 1 in which the sites 
ki, . . . ,kp are occupied. In effect we have removed from the lattice the site 
lying to the right of each occupied site. Note that our earlier inequality 
2p < n becomes, in terms of n', simply p < n' — 1, which reflects the fact 
that the new lattice has only n' — 1 sites. The concurrence is still given by 
C = (2/n)|y|, where 

P r -1 

y = J2 J2 J2 Ki,...,kph[,...,k'Jk'^,k,+i n ^K,kr ■ (24) 



q=l ki<-<kp k[<---<k'p 



We can express y more simply by introducing creation and annihilation 
operators for each site. We associate with site k the operators 



Cfc = „ „ and 4 = K „ , (25) 





which are represented here in the basis {|0), |1)}. In terms of these operators. 
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we can write y as 

n'-2 

y = (CIE4c.+ilC)- (26) 

fc=l 

Our problem is beginning to resemble the nearest-neighbor tight-binding 
model for electrons in a one-dimensional lattice. The Hamiltonian for the 
latter problem — assuming that the spins of the electrons are all in the same 
state and can therefore be ignored — can be written as^ 

n'-2 

^ = - E(ck+i + 4+iCfe), (27) 
fe=i 

where we have taken the lattice length to be the same as in our problem, 
namely, n' — 1. From Eqs. (26) and (27) we see that (C|i?|C) — — 2Re(y). 
This expectation value is not quite what we need for the concurrence: the 
concurrence is proportional to the absolute value of not its real part. 
However, as in our earlier example, for the purpose of maximizing C there 
is no advantage in straying from real, non-negative values of bki,...,kp- If we 
restrict our attention to such values, then the absolute value of y is the same 
as its real part, and we can write the concurrence as 

C=~{C\H\C). (28) 

Thus, maximizing the concurrence amounts to minimizing the expectation 
value of H, that is, finding the ground state energy of the tight-binding 
model, as long as the ground state involves only real and non-negative values 
of hi,...,kp- 

The onc-dimcnsional tight-binding model is in fact easy to solve [16, 17]. 
Its ground state is the discrete analogue of the ground state of a collection of 
p non-interacting fermions in a one-dimensional box. In our case the "walls" 
of the box, where the wavefunction goes to zero, are at A; = and k — n\ 
and the ground state |Co) is given by the following antisymmetrized product 
of sine waves: 



hk, fe„ oc A 



sm 



^In Eq. (27) the operators c and are fermionic, whereas those defined in Eq. (25) are 
not, because they do not anticommute when they are associated with different sites. We 
could, however, use our c's to define genuinely fermionic operators in terms of which the 
extremization problem has exactly the same form [16]. 
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Here A indicates the operation of antisymmetrizing over the indices ki, . . . ,kp. 
In the range of values we are allowing for these indices, that is, < fci < 
k2 < ■ ■ ■ < kp < n' , the coefficients bki^,.,,kp are indeed non-negative, so that 
Eq. (28) is valid. 

The ground state energy, from which we can find the concurrence, is 
simply the sum of the first p single-particle eigenvalues of H. There are 
exactly n' — 1 such eigenvalues, one for each dimension of the single-particle 
subspace; they are given by 

= -2cos(^), m= 1. (30) 

Thus the concurrence is 

C^-i(Co|H|C.>^^Eco.Q. (31) 

Doing the sum is straightforward, with the following result: 

cos(p7r/n') — cos((p -|- 1)71/11') + cos(7r/n') — 1 

1 - cos(7r/n') J' ^ ' 

Recall that n' — n — p + 1. Eq. (32) gives the largest value of C consistent 

with our conditions, for fixed values of n and p. Note, for example, that 
when n = 5 and p = 2, Eq. (32) gives C = V2/5, just as we found before for 
this case. 

We still need to optimize over n and p. It is best to make the block size n 
very large — any state w that is possible with block size n is also allowed by 

block size 2n — so we take the limit as n goes to infinity. Let a be the density 
of occupied sites — that is, a = p/n — and let n approach infinity with a held 
fixed. In this limit, the concurrence becomes 

Ciim = -(!-«) sin (-^). (33) 

Taking the derivative, one finds that Cum is maximized when 

tanf- ) = - , 34 

which happens at a = 0.300844, where Cum = 0.434467. This is the highest 
value of concurrence that is consistent with our method of constructing the 
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state of the chain and with our three conditions on |,^). Note that it is 
considerably larger than what we got in our first example, in which a string 
of singlets was mixed with a shifted version of the same string — one might 
call this earlier construction the "bicycle chain" state. Unlike the bicycle 
chain state, our best state breaks the symmetry between the basis states |0) 
and |1): the fraction of qubits in the state |1) is about 30% rather than 50%. 
Of course the entanglement would be just as large if the roles of |1) and |0) 
were reversed. 

It is interesting to ask what value of entanglement of formation the above 
value of concurrence corresponds to. As a function of the concurrence, the 
entanglement of formation is given by 



where h is the binary entropy function h{x) = — [a;log2a; + (l— x) log2(l — x)]. 
For the above value of concurrence, one finds that the entanglement of forma- 
tion is Ef = 0.284934 ebits. (For the bicycle chain state, the entanglement 
of formation between neighboring pairs is only 0.118 ebits.) 

If one can prove that this value is optimal, then it can serve as a reference 
point for interpreting entanglement values obtained for real physical systems. 
A string of spin- 1/2 particles interacting via the antiferromagnetic Heiscnberg 
interaction, for example, has eigenstates that typically have some non-zero 
nearest-neighbor entanglement. It would be interesting to find out how the 
entanglements appearing in these states compare to the maximum possible 
entanglement for a string of qubits.^ 

Clearly the problem we have analyzed here can be generalized. One can 
consider a two or three-dimensional lattice of qubits and ask how entangled 
the neighboring qubits can be. If we were to analyze these cases using as- 
sumptions similar to those we have made in the one-dimensional case, we 
would again find the problem reducing to a many-body problem, but with 
less tractable interactions. Assuming that pairwise entanglement tends to 
diminish as the total entanglement is shared among more particles, one ex- 
pects the optimal values of C and Ef to shrink as the dimension of the lattice 
increases. 

*Since the original version of this paper was written, the question about the antifer- 
romagnetic Heiscnberg chain has been answered for the ground state [18]: though the 
nearest- neighbor concurrence of the ground state is high (C = 0.386), it is not optimal. 
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